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EARLY TIME COUPLING STUDIES USING A ID HYBRID CODE 

I. INTRODUCTION 

The understanding of atmospheric effects generated by a HANE is of 

vital importance to DNA. It is well known that HANE's cause severe 

3 

disturbances in the atmosphere and ionosphere which impact C I systems. In 
particular, long lasting, large scale ionization irregularities degrade 
radar and communications signals. To develop predictive capability for 
such HANE effects, DNA sponsored a major research effort at NRL leading to 
theoretical and computational models of HANE. Significant progress was 
made in this area during the past decade leading to the development of 
sophisticated raultifluid codes which incorporate in a self consistent 
manner anomalous transport due to instabilities. Recent advances in both 
plasma theory and numerical techniques allow us today to produce 
substantially more refined and accurate models of the early time HANE 
phenomenology. These advances indicate the need for a better incorporation 
of kinetic effects due to strong deviations of the particle distributions 
from Maxwellian. The new codes allow us to compute the laminar and 
turbulent dynamic ion response over long time and space scales and for a 
mixture of ionized species. In this way laminar, electrostatic, Larraor and 
turbulent coupling processes can be followed dynamically and their relative 
strength and dynamic interplay can be assessed. This paper is the first to 
describe this series of investigations, focuses on a particular process, 
early time debris-air coupling, and stresses quantitative understanding of 
the underlying physics. 

Early time coupling had been a controversial subject in the DNA HANE 
community in the 1970's. Two "opposing" views were maintained: NRL, LASL 
and ARA advocated short-range coupling associated with plasma turbulence, 
while MRC argued for "Larmor coupling" associated with the ion gyration in 
the ambient magnetic field. To a certain extent, the issue was not whether 
or not these coupling mechanisms worked, but which mechanism was 
dominant. The turbulent coupling process developed at NRL was studied by 
incorporating theoretical estimates of anomalous transport coefficients in 
a ID multi-fluid code. Although this procedure yielded results consistent 
Manuscript approved February 14,1985. 
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with HANE data and provided strong support for turbulent coupling, it was 
not a self-consistent technique. Moreover, it was not capable of handling 
kinetic effects such as caused by reflected ions; such reflected ions are 
known to occur in high Mach number shocks. The present study, using a ID 
hybrid code, overcomes these deficiencies. Namely, it allows instabilities 
to be excited which affect the plasma evolution in a self-consistent 
manner, i.e., no model anomalous transport coefficients are used, and 
reflected ions are permitted. Furthermore, the ID hybrid code inherently 
includes "Larraor coupling" so that a direct comparison of the relevant 
importance and interplay of the coupling mechanisms can be made. 

The organization of the paper is as follows. In the next section we 
present an overview of the basic physical processes occurring as debris 
streams through the air. In Section III we present results of two 
simulations which demonstrate the various aspects of coupling. Finally, in 
the last section we summarize our results and discuss their implications to 
HANE phenomenology. 

II. OVERVIEW OF PHYSICAL PROCESSES 

The evolution of the debris-air plasma is characterized by three 
temporal stages: (1) magnetic field compression; (2) piston formation; and 
(3) shock formation and evolution. The bulk of this paper describes in 
detail piston dynamics; for completeness we briefly discuss the other two 
stages. 

In the initial stage, the debris streams outward from the burst point 
and picks up air electrons which merge with the debris electrons to form a 
single electron distribution. The compression of the air electrons leads 
to a compression of the magnetic field since the air electrons are "tied" 
to the field lines due to the low value of resistivity. This is the first 
stage of the debris-air interaction and has been well-documented 
theoretically (Longraire, 1963; Sloan, 1970; Lampe and Hernandez, 1972), 
numerically (Clark et al., 1974), and experimentally (Ripin et al., 1984). 

The next stage involves the dynamics of piston formation and will be 
discussed in detail in the following section. Brieflv, this stage is 
dominated by the interaction of the two ion streams, i.e., the debris ions 
and the air ions. It is this phase where various coupling mechanisms 
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become apparent. The debris ions slow down, while the air ions are picked 
up; there is momentum exchange between the ion species. There is also a 
conversion of the directed kinetic energy of the debris to thermal energy 
of both ions and electrons, i.e., heating and acceleration occurs. And 
finally, an electrostatic potential forms at the leading edge of the 
magnetic field compression which causes some debris ions to be accelerated 
ahead of the magnetic piston and some air ions to be reflected. This final 
point, although well-known, can be very important to the formation of 
debris patches in the conjugate regions since it results in very high 
velocity ions. 

Lastly, the final stage involves the formation of a shock wave which 
"runs" ahead of the piston. We do not describe this phase in detail in 
this paper, but do show that a minimum time is required for high Mach 
number shock formation. This minimum time is basically the time needed for 
the shock to "run" ahead of the accelerated debris or reflected air ions. 


III. SIMULATIONS AND RESULTS 


Prior to discussing the details of the simulations, we present a 
simple model equation that elucidates the various coupling mechanisms to be 
studied. Moreover, it highlights the distinction between previous models 
(e.g., KLYSMA) and the present work. We consider the following momentum 
equation in the x direction for an air ion (i.e., radial direction 
perpendicular to the ambient magnetic field B_■ Be^). 


dV 


eE. 


dt 


~ - v.) 


a ay 


a) 


where the subscript a refers to air, B/m^c, v* is an anomalous ion- 
ion collision frequency, = ^d e x stream * n ® debris velocity, e is 
the charge, m is the mass, and Z is the charge state. The first term on 
the RHS of (1) arises from a laminar electric field usually associated with 
the Leading edge of the magnetic compression. It acts to accelerate debris 
ions and to reflect air ions. The second term in (1) is the magnetic force 
which is associated with Larmor coupling. The final term in (1) 
corresponds to turbulent "pick up" of the air and arises because of plasma 
instabilities. For plasma turbulence such that v* > it is clear that 
turbulent coupling can dominate over Larmor coupling. 
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Notice that the force ^ a V associated with Larmor coupling is 

proportional to the value of V a ^ which, before any form of coupling occurs, 

is V = V ■ 0. Therefore the laminar force eE /m and the turbulent 

*** x a 

force v *|V ax - V^j ~ v*V^ will dominate initially even for v* « & a (i.e., 
as long as v* > J2 a (V a y/V^)). The time evolution of V a ^ will be given by 


In deriving eq. (la) we assumed an infinitely conducting (n -*• 0), low 8 
(8 << l) situation so that E = V^B/c. From eqs. (1) and (2 ) it is 
obvious that the time scale for significant Larmor coupling is R & t > 1, 
i.e., the value of V a ^ should approach V^. At earlier times the laminar 
and turbulent forces will dominate. The turbulent forces will also 
dominate even at t » & a ^ as long v* > $ a jV ax - V^J/V^. Another important 
point follows from (1) and (2) if we neglect the laminar and turbulent 
forces. If L is the size of magnetic field compression then as long as V ax 
< Vj the maximum value of V a y will be = Cl a L. As ment l° nef * previously 

the Larmor momentum coupling force in the x direction becomes significant 
if V ay 3 fl a L. Therefore for effective Larmor coupling L > V d /fl a . Namely, 
an extremely broad magnetic field profile will be required. We will 
further comment on this point later on. 

In previous work, theoretical estimates of v* were incorporated in a 
multi-fluid code (KLYSMA) which solved momentum equations similar to (1) 
and (2) (Clark et al., 1974). However, the turbulent interaction involves 
ions interacting with fluctuating electric fields, i.e., v* 3 \>*(5E). The 
present simulations incorporate such effects self-consistently. The air 
momentum equation solved, analogous to (1), is 


+ SI V 
a ay 


where is self consistently generated because of plasma instabilities 
(e.g., ion-ion streaming instability (Papadopoulos et al., 1971) and is 
allowed to act directly on the ion motion, rather than through a model 
equation such as (l). Given this background, we now present the simulation 





















The simulations were performed using a one dimensional quasi-neutral 
hybrid code (Chodura, 1975; Winske and Leroy, 1984) whose description is 
given in Appendix I. The code incorporates ion kinetic effects by 
following the trajectories of "superparticles", numerically representing 
many actual ions, in 3 dimensions in velocity space and one dimension in 
configuration space. The electrons are approximated as a massless fluid 
described by momentum and energy transport equations; the electron pressure 
is assumed to be isotropic. These equations, with Maxwell's equations, are 
solved on a one dimensional spatial grid using moments calculated from the 
ion distributions. Anomalous resistivity, resulting in magnetic field 
diffusion and electron ohmic heating, is included as a parameter in the 
equations with a value consistent with that expected from current driven 
cross-field instabilities. 

The boundary conditions in the code correspond to the conditions we 
expect at early times. A dense ion beam is injected continuously into the 
simulation box which is initially filled with a tenuous stationary 
plasma. In our present results both ion populations are H + , but other ion 
species can be easily considered and will be presented elsewhere. We set 
up a finite width magnetic field compression within the leading edge of the 
ion beam. This represents the field compression expected from the cross¬ 
field beam motion which, because of finite resistivity, will diffuse into 
the beam. The width of the compression has been varied to study various 
physical phenomena. The magnetic field is held fixed at the edges of the 
simulation box. 

We present below the results of two simulations selected to illustrate 
the underlying physical coupling processes. The initial state is typical 
of the piston formation stage. It is shown in Fig. la, for the case of M A 
= V/V a * 4, where is defined on the basis of the ambient (i.e., 

upstream) parameters. A magnetic pulse consistent with B/n = constant is 
introduced at the left hand boundary. In this paper we do not study the 
formation of the pulse but assume its nature in accordance with previous 
work mentioned in Section II, and its width L to illustrate the important 
physics. In the figures the velocity is in terms of the upstream value of 
V A and of time in terms of the upstream value of £? ^. The units of length 
are in terms of the ion cyclotron value computed for air with velocity 


but with the value of the downstream (i.e., compressed) magnetic field. 
This representation has been selected for clarity of the underlying 
physics. We consider (1) a "broad" Lg = 2R^ (R^ in upstream values) pulse 
with ^ 1 4 and n a /n a = 16 and (2) a "narrow" Lg = pulse with * 8 

and n 0 /n a * 64. Here V^, n^ are the debris velocity and density, n a the 
air density and is defined on the basis of the upstream conditions. 

A. Broad Pulse (Lg * 2r^) 

This simulation is initialized with a magnetic pulse width L B ~ 2r Li 
V^/V^ = 4 and n ( j/n a = 16. We show a series of results for this run in Fig. 
1. Figure 1 displays four important quantities versus distance. The 
bottom curve represents the magnetic field magnitude, the middle curve (or 
"dots") represent particle velocities of the debris and air ions (the 
debris ions have ~ 4V^ while the air ions have ~ 0), and the top 
curve represents the electric field in the x direction. Figure la is at t 
= 0.0059 where is based on the upstream value of B. Note the broad 
pulse in B with sharp gradients at the leading and trailing edges. 
Associated with these gradients are strong laminar electric fields at these 
edges as shown in the top portion of Fig. la. These fields strongly affect 

the ions. Note in the middle panel that some debris ions are being 

accelerated ahead of the B pulse, and air ions are being picked up in a 
reflection process. 

Figure lb depicts the system at t = 0.0386 and several interesting 
features are developing. First, the B field pulse has broadened 

considerably. The leading edge of the pulse has been "dragged" out to 

x ~ 4r^, and the gradients are not as sharp. Second, the laminar electric 
field has fallen in intensity and is accompanied by a considerable amount 
of turbulence in the region 0 < r/r^ < 2.4 caused by of the ion-ion 
streaming. This turbulence causes some momentum exchange between debris 
and air ions, and thermalization of the particles in this region. Finally, 
also at this time, a fraction of debris ions have been accelerated to a 
velocity greater than twice the initial velocity. 

Figure lc shows the system at t = 0.0874 0^ . The magnetic pulse 

continues to broaden and now extends to x ~ 6r Li . A laminar electric field 
still exists at x ~ 3.6 However, the turbulence in the region 0 < 





x/r^^ < 3.6 has subsided since the relative velocity between debris and air 
ions is small. The ion motion in this region is dominated by the magnetic 
field and, in effect, is where "Larmor coupling" is taking place. Perhaps 
the most interesting region is 3.6 < x/r^ < 6.0. It is clear that debris 
acceleration and air reflection are continuing. Furthermore, turbulence is 
developing at the debris-air interface (x ~ 4.0 r-^). 

Finally, Fig. Id shows the system at t = 0. 1542 £5^*. The magnetic 

pulse is quite broad now, extending out to x ~ 8.0 r^. Very intense 

electrostatic turbulence is occurring at the debris-air mix in the region 
4.2 < x/r^ < 8.2 which is producing coupling between the ion species and 
therraalizing the particles. Debris particles have been accelerated up to 

V ~ 3V^. In the region 2.0 < x/r^ £ 4.0 tbe a ^ r i ons have a velocity 

V ~ 5.0 V^ while the debris ions have a velocity V ~ 3.0 V^. In this 

region the momentum coupling is through the magnetic field (Larmor 

coupling) and the debris and air ions perform gyrations about the magnetic 

field. Very little turbulence is evident in the electric field, confirming 

the fact that v* < fl . 

a 

Thus, in this simulation it is found that three important processes 
are affecting coupling, energization, and thermaiization: (1) laminar 
electric fields; (2) turbulent electric fields; and (3) magnetic field 

pick-up. The laminar fields are most prominent early In the run (Fig. la) 
and cause debris acceleration and air reflection. This is crucial to early 
time HANE processes since debris ions can be accelerated up to ~ 3 times 
their initial expansion velocity, and air ions are "picked-up". Turbulent 
electric fields cause coupling and thermaiization of the ion species. In 

this run they are most effective at the debris-air interface, and ahead of 
the interface (Figs. 1c and Id). Finally, "Larmor coupling" is observed 
within a section of the broad magnetic pulse in Figs, lc and Id, i.e., the 
debris and air ions have rotated about the B field so that the air ions 

have a larger v velocity than that of the debris (V flX > Vj x ). This process 

supplements any incomplete part of the turbulent coupling and operates on a 

long time scale. 
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B. Marrow pulse (Lg - 0.1 r^) 

This simulation is initialized with a magnetic pulse width -0.1 
r^j and V^/V^ = 8 and n ( j/n a = 64. The results are shown in Fig. 2. In 
Fig. 2a we show the plasma system at t = 0.0206 . Already by this time 

a substantial number of air ions are being reflected from the laminar 
electric field at x ~ 1.5 r^j* Furthermore, debris ions have been 

accelerated up to V ~ 2.5 and extend out to x ~ 4 r^. In this region 
there is a considerable amount of electrostatic turbulence as seen in the 
electric field plot. Figure 2b is at t = 0.0411 . The magnitude of the 

turbulence electric fields have increased from the previous time shown and 
are quite intense at the debris-air interface (x ~ 3.0 r^). The air ions 
are being "picked up" because of their turbulence and thermalized. The 
magnetic field has been dragged out to x ~ 7.0 r^ by accelerated debris 
ions. Finally, Figs. 2c and 2d show the ions and fields at t = 

0.0617 and t = 0.0771 respectively. The general features are 

similar to t = 0.0411 (Fig. 2b). There are intense turbulent fields at the 
debris-air interface which transfer momentum from the debris to air ions, 
and also thermalize the particles. Debris ions extend out to x > 10.0 r^ 
by t = 0.0617 as well as the B field and low level turbulence. In 

this run there is no indication of Larmor coupling. One interesting 

feature seen in Fig. 2d is the "structuring" of the magnetic pulse. In 
particular, there Is a sharp gradient that has developed at x ~ 4.8 r-^ 
which produces a laminar electric field. This field acts to energize air 
ions to such that V ~ 1.8 at x ~ 5.0 r-^. 

IV. DISCUSSION 

We have presented the results of simulations of very early time 
processes (i.e., t < 12^) using a ID hybrid code. The purpose of this work 
has been to elucidate important physical phenomena rather than to perform 
an actual early time simulation of a HANF. The focus of this studv has 
been primarily on coupling mechanisms; an area that has been somewhat 
controversial within the DNA community. We have found that there is no 
"single" coupling mechanism but that several mechanisms can be effective at 
various times throughout the debris expansion phase. In the first 
simulation, initialized with a "broad" magnetic pulse, the initial coupling 











occurs through laminar electric field at the leading edge of the pulse. 
Subsequent to this, turbulent coupling occurs at the initial debris-air 
interface while Larmor coupling occurs within the main body of the magnetic 
pulse. In the second simulation, initialized with a "narrow" magnetic 
pulse, the only coupling observed was because of the electric field 
(laminar and turbulent). 

A very important result of this study is the generation of very 
energetic debris ions. This is caused by the laminar electric field at the 
leading edge of the magnetic pulse. A fraction of the debris ions can be 
accelerated up to roughly three times the initial debris expansion velocity 
(i.e., ~ 3Vj). It is well-known that the initial formation of ion 
debris patches occurs on a time scale faster than that calculated from the 
initial expansion velocity. Thus, the accelerated debris ions observed in 
the simulations is consistent with HANE observations. 

Finally, this study dramatically highlights the importance of kinetic 
phenomena for the understanding of early time processes. Kinetic effects 
such as plasma turbulence, reflected and accelerated ions, and 
therraalization of the ions are not fluid processes, and cannot be recovered 
using only MHD codes. Detailed studies are necessary using codes such as 
the one described in this paper to properly understand the dynamics of 
early time HANE expansions. 
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(a) 

Fig. 1 — Evolution of the ambient magnetic field (bottom curve), debris and air ions (middle 
curves/dots), and electric field in x-direction (top curve) as a function of time. The parameters are 
M a - VJV A - 4, L b ~ 2r u , and njn a - 16. (a) t = 0.0059 ftf 1 ; (b) t - 0.0386 n,"*; (c) t * 
0.0874 n: 1 ; and (d) t - 0.1342 nr 1 - 
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Fig. 1 (Corn’d) — Evolution of the ambient magnetic field (bottom curve), debris and air ions 
(middle curves/dots), and electric field in x-direction (top curve) as a function of time. The 
parameters are M A - V d / V A =» 4, L a — lr Ln and rtjn a - 16. (a) t = 0.0059 fi~'; (b) t = 0.0386 
nr 1 ; (c) t * 0.0874 nr 1 ; and (d) t = 0.1342 fl 
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(c) 

Fig. I (Corn’d) — Evolution of the ambient magnetic field (bot.om curve), debris and air ions 
(middle curves/dots), and electric field in x-direction (top curve) as a function of time. The 
parameters are = VJ V A - 4, L g — 2 r Li , and njn a = 16. (a) t = 0.0059 ft (b) t - 0.0386 
fir 1 : (c) t = 0.0874 nr 1 : and (d) t - 0.1342 ftr 1 . 
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Fig. 1 (Cont’d) — Evolution of the ambient magnetic field (bottom curve), debris and air ions 
(middle curves/dots), and electric field in x-direction (top curve) as a function of time. The 
parameters are M 4 * VJ V A = 4, L g — 2 r Li , and njn a - 16. (a) t = 0.0059 (b) t * 0.0386 

nr 1 ; (0 t - 0.0874 nr 1 ; and (d) t = 0.1342 nr 1 . 
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APPENDIX 

HYBRID CODE DESCRIPTION 


We describe here the numerical details of the hybrid code used in this 
paper. This description is taken from Winske and Leroy (1983) with minor 
changes necessary due to the HANE application. Some of the techniques 
(e.g., particle-in-cell method) are standard and well known, while others 
(e.g., solution of the field equations) have been applied primarily to 
fusion-related problems and hence are generally not familiar to the space 
science community. The discussion is divided into five parts: (A) an 
overall description of the simulation methods, (B) treatment of the ions, 
(C) solutions of the field and electron fluid equations, (D) refinements to 
the basic model, and (E) discussion of the initial and boundary 
conditions. The methods pertaining to the first four parts are rather 
general and could thus be applied to any problem where a 'hybrid' (kinetic 
ion, fluid electron) description is appropriate. Specific adaptation of 
the method to the HANE study enters primarily through the last topic, the 
initial and boundary conditions. 

A. Overall Description 

The simulation method is particle-in-cell (Morse, 1970) with the 
electrons treated as a fluid (i.e., cell) quantity rather than as discrete 
particles. The assumption of a massless, charge neutral electron fluid 
thus eliminates the restriction of time and spatial lengths to he inverse 
electron plasma frequency and electron Debye length, usually associated 
with full particle codes. The hybrid models discussed here were originally 
developed for and extensively applied to dynamical studies of high density 
pinch experiments by Chodura (1975), Sgro and Nielson (1976) and Hamasaki 
et al. (1977). These models apply strictly to the case of only one spatial 
dimension; extension to two dimensions is possible, although nontrivial, as 
discussed later. Even with only one spatial dimension (which is along the 
x-axis and thus implies 9/9y = 9/9z = 0), all three velocity, magnetic 
field and electric field components are included. The simulation region 
has length L, divided into N cells, each of length Ax = L/N. There is an 
additional ('ghost') cell at each end of the system, which helps in keeping 
track of the particles entering and exiting the system and in setting up 


the boundary conditions. The cell quantities (i.e., fields, electron fluid 
properties, ion velocity moments) are specified at the cell centers. 
Because the restriction to electron spatial and temporal scales has been 
eliminated, much larger time steps can be used. Typically, the time 
step (At) is limited by the condition that ion gyromotion is well defined 
and the cell size is constrained by the condition that the fastest ions do 
not traverse one cell in one time step. The cell size and system length 
are also chosen to resolve length scales of interest to the problem. 

3. Dynamics 

The ion component is modeled by a discrete set of particles. The ion 
distribution is advanced in time by stepping forward each particle in time 
under the influence of the local, self-consistent Lorentz force. The 
motion of the ions in a four-dimensional phase space (v x , v , v z , x) is 
solved by the particle-in-cell technique, using a second-order-accurate but 
non-reversible scheme (Nielson and Lewis, 1976). The equations for the ion 
advance (with superscripts denoting the time level) are: 


—0 

v x B 


y 3 = y' 1/2 + (h/2) E° 

v 1/2 - fv _l/2 + h(E° + gB° + + P°) 

1/2 


X 3 = X 3 + it V 


(Al) 


where h = At e/m , f = l-fh 2 /2)B° • B°, g = fh/2)(v~ 1/2 • B°), e is the 

charge and m^ is the mass of the ions. The electric [EJ and magnetic 

fields (b) are evaluated at the particle position x^. P is a mean friction 

force, P = -eg • J, exerted by the electrons as a macroscopic force 

only, J is the current and g represents a phenomenological anomalous 

resistivity which gives rise to Ohmic heating, as will be explained 

later. After all of the ions have been advanced, the ion density fn^) is 

updated, by averaging over the positions of all the particles in each 

cell. The ion mean velocities fv. 1 are also needed at time level 1. The 

^~i - 

velocities are pushed ahead one half time step, 

-1 1/2 1 ~ l/2 X ~ l 1 

v 1 =■ v 1/2 + (h/2)(E 1 + ---+ P 1 ), (A2) 
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using the fields at the new particle positions x\ in order to accumulate 
the velocity moments V . 

The assumption of charge neutrality then implies 


and the continuity equation gives 


jL 

3t 


(n + n, ) = 0 
^ e i 1 


3x 


3 _ 

3x 


en (V 


ix 


V 1 


ex' 


(A3) 


(AA) 


or 


V 

ex 


ix‘ 


(A5) 


Thus, the electron density and one component of the electron fluid velocity 
are determined. 


C. Field-Fluid Equations 

With the assumption of zero electron mass the electron momentum 
equation reduces to 


dV 


n ra 
e 


~e 
e dt 


0 


(E 


V 

~e 


x B 


: ) - 


3_ 

3x 


P * 


n P 
e~ 


( A6 ) 


where p g is the scalar pressure and is the velocity of the electron 
fluid. The transverse (y and z) components of (A6) can be written as Ohm's 
law 

, V x 3 

J = g • (E + —— ). (A7) 


In general, the magnetic field is inclined at some angle with respect to 

the x axis: 3 = |b| (b , b , b ) with b^ + + b^ = 1, and if a. fa.) 

~ x y z' x y z 111 

represents the conductivity parallel (perpendicular) to the magnetic field, 
the conductivity tensor is (Krall and Trivelpiece, 1973): 
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2 s 3 


-l 


a. a a 

1 y z 

a o~ a , 

y 2 x 

a a a. 

zx3 


(A8) 


where 


a, = a. + fa - a. lb b 

1 1 fl 1- x x 

a 0 = a. + fa - a.lb b 

2 1 i) 1- y y 

a- = a, + fa - a.lb b 

3 1 v H 1 ; 2 2 

a = fa- a, lb b 
x '■I! y z 


,a. - a.Jb b 
1 n L J x 2 


a = fa. - a. lb b . 

2 H 1' x y 

The transverse fields can be expressed in terms of the vector 
potentials in the usual manner fB = * A; E * - fl/c)f3A/3t)1: 

3A 

\ ’ 3T 

3A 

B z ’ 3 ^ 


T 3A 

I _I 

c 3t 
. 3A 

_L __jl 

c 3t 


(A9) 


and 7*3=0 implies 


B = constant. 


(A10) 


(The last field component, E x , is derived from the x-component of (A6) as 
will be shown Later.) 
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Neglecting the displacement current, Ampere's law can be expressed 


3“ A 


4 TT 


3x 


9 

3“A 


(All) 


3x 


Substituting in Ohm's law (A7) and the definitions of the field 
components, Eq. (11) becomes a coupled system of equations: 


3~A 

y 

= f 

3A 

f -JL 

3t » 

3A 

z 

3A 

3A 

3x 2 

r i 

3t ’ 

’ 3x ’ 

1 3x j 

9 

3“A 

z 

= f 

3A 

f— l 

3t ’ 

3A 

z 

3A 

_ L 

3A 

—L ), 

' 3x - 

Cl 2 

r 2 

3t ’ 

1 3x ’ 

to 

solve 

this 

system a 

fully 


(A12 ) 


difference scheme is used (Richtmyer and Morton, 1967). If we let u. = A^ 


at the jth cell position and n-th time level (and w 


A z in the same way). 


the left hand side of (A12) is differenced as follows: 


3 2 A 


n+L , n+1 _ n+1 

u... + u. . - 2u, 

3+1 3-1 3 


(A13) 


3x 


Ax 


while on the right hand side 

n+1 n 


3A 


y _ 


3t 


3 


At 


and 



n+1 n+1 

u.,, - u . , 
ljl 3~t 

2 Ax 


(A14) 


When this is substituted into (A12), the resulting system of equations can 
be written symboLicaLly as 



















X. 

~3 


n+l 
a. 

J 

n+l 

w. 


(A16) 


and A^, B^, and Cj are 2x2 matrices that depend only on components of a and 

V (evaluated at cell j at time level n) and constants, while D also 
~e ~ 

depends on u.., w^. Thus, all the coefficients A, B, C and D are explicitly 
known. Note that there are N cells in the computation mesh j**l, 2, ... N 

with ghost cells (.1*0 and i=«N+l) at each end. 


The set of equations (A15) is easily solved. Assuming a solution of 
the form 


X. = E. • X.^. + F. j = 0, 1, ... N (A17) 

and substituting it into (A15), it then follows that 

E. ■ - Tb. + C, • E. . I' 1 A. 

“J L *J a 1 a ]-l J 3 3 

Lj - [Bj + C. • E..J' 1 • [D. - C i • F.^1 (A18) 

Because all the A., B., C. and D. are known, if E„ and F~ are known (from 

a 3 *3 3 j ~3 ra 0 ~0 

boundary conditions, as explained Later), (A18) can be solved in ascending 

order to obtain E., F. ; E„, F 0 , .... E„, F, T . Then if X ,, can be 
«1 ~1 «2 ~2 »N ~N ~n+l 

determined (again from boundary conditions), (A17) can be solved in 

descending order to obtain X^, X^j, • ••, Xq* 

Once the components of the vector potential are known, the new values 
of the electric and magnetic field follow from (A9) and the components of 
the current from (AIL), again using the differencing in (A13) and (A14). 
The transverse components of are then easily found from the current and 
the ion velocity moments: 

V . = V., - J./en } = y or z. (A19) 

«3 ij 3 
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The last electron fluid quantity to be calculated is the temperature. 
The electron energy equation can be written as 


2 3_ 
2 3t 


(n T ) 


_a_ 

3x 


[=r n T V ) + 
2 e e ex' 


3V 

_ ex 
n T —,— * Q. 
e e 3x 


(A20) 


2 

The source term Q includes resistive heating (nJ ) and loss mechanisms, 
such as thermal conduction or radiation, depending on the application. 
Letting 


n T , 


(A21) 


(A20) can be rewritten in the form (with y » 5/3) 
. . 3V 

ri_ + \j 

'■at 


ex 3x J e 


- YP 


ex 


e 3x 


+ (y-1)Q 


(A22) 


The differencing for 3p /3t and 3V /3x is identical to (A14); for 

e ex 

stability reasons "donor cell" differencing is used for the convective terra 
(Richtmyer and Morton, 1967): 


, Hi 

ex 3x 


exj 

Ax 


n+1 

n+1 > 

V , > 0 
exj 

ej 

P e.i-2^ 

n+1 

n+1 1 

V . < 0 
exj 

e.i+1 

- V) > 


(A23) 


Again, a tridiagonal system is obtained. 


a jVi + B j x j + “J 




(A24) 


except that now ail quantities are simple scalars. As before, an assumed 
solution of the form 


leads to 


x j - e jVi + F i 


tc iVi r S 


(A25) 
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(A26) 


Tb, 


1 tc iVi r ''"i - c j p j-i-* 


Appropriate boundary conditions lead to a determination of Eq, F q (and then 
to all Ej, Fj in ascending order, using (A26)) and X^ + ^ (and thus to all Xj 
in descending order through (A25)). Since n ( a n^) is known, T e is then 
obtained from (A21). An alternative to (A22) which is sometimes useful 
(Sgro, 1978) is to write the energy equation (A20) in terms of the 
entropy [s = in(p e n e ^)1 instead of the pressure. 

Finally, the x component of the electron momentum equation (A6) can be 
solved for E x : 


E 


x 


c 




x B 
~x 


) ~ — 
' n o 



(A27) 


This electric field is needed to maintain charge neutrality. 

The computation loop through one time step can thus be summarized as 
follows: 

1. Advance the ions one time step and then calculate the ion velocity 
moments (n^ ■ n g , V ix » V ex. V iy. V i«>- 

2. Compute the plasma conductivity g, according to some prescription. 
(The resistivity rj is g ). 

3. Solve the coupled equations for A„, A_; we can then easily 

y c 

calculate E y , E z , B y , B z (B x = constant). From the currents 
(Ampere's law) and the ion moments calculate the other two 
electron velocity components, V gy and V gz . 

4. The electron temperature is obtained next from the solution of the 
differenced energy equation. 

5. E x is then calculated from the x-component of the electron 
momentum eauation. 

6 . Thus, all field and electron fluid components are known and we are 
ready to move the ions again. 
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D. Refinements 


The simulation scheme described thus far is rather general, the only 
assumptions being quasineutrality, zero electron mass and one spatial 
dimension. All of these conditions can be replaced by more appropriate 
ones, as the physical situation dictates. For example, for the study of 
low frequency ion waves an adiabatic electron model is more appropriate; 
such a model has been successfully used by Okuda et al. (1978). In some 
situations phenomena in the lower hybrid frequency range are of interest. 
In this case electron inertia effects are non-negligible; they have been 
included in the one-dimensional model of Liewer (A1976) and the two- 
dimensional model of Hewett and Nielson (A1978). Two-dimensional 
simulation models with m g = 0 have been successfully employed by Byers et 
al. (1978), Hewett (1980) and Harned (1982). 

The HANE simulations described in this paper take the conductivity to 
be constant. A further level of sophistication is to include anomalous 
processes due to the microphysics, which is occurring on time and distance 
scales shorter than those resolved in the hybrid model, by means of more 
complicated transport coefficients. This is done by expressing the 
conductivity (or resistivity) as a sum of two terms, one representing 
classical (Spitzer) effects and the other due to anomalous effects arising 
from microinstabilities due to cross-field currents. Two types of 
anomalous transport coefficients have been used successfully to model the 
behavior of laboratory plasmas: one type uses a semi-empirical expression 
(Chodura, 1975; Sgro and Nielson, 1976), the other is based on a 
quasilinear analysis of known instabilities (Davidson and Krall, 1977). 
The transport processes can be further refined to include electron thermal 
conductivity (Sgro, 1980) and anomalous ion heating (Hamasaki et al., 
1977). Multispecies ions can also be included in a straightforward manner 
(Sgro, 1980; Sgro and Winske, 1981). 






E. Initial and Boundary Conditions 

In order to simulate the piston formation stage of the HANE problem, 
appropriate initial and boundary conditions have been implemented in the 
hybrid code. We assume initially that, a dense cool debris plasma flows 
across the magnetic field into a tenuous stationary air plasma. In this 
paper, we assume the debris and air to consist of the same ion species. 
This assumption is not necessary and we are presently investigating the 
effects of different debris composition. We assume further that there is a 
magnetic compression of finite width imbedded in the leading edge of the 
debris stream, formed as a result of the preceding magnetic compression 
stage of the coupling is simply maintained at the right boundary. Ions of 
all species are effectively absorbed at the boundaries. The final boundary 
condition is on the magnetic field which is held constant at the ambient 
(air) value at both boundaries. The magnetic compression is taken to be 
proportional to the electron compression. The width of the magnetic 
compression is an input parameter in the code. Thus initially in the 
simulation space we establish a debris ion stream extending a finite 
distance from the left boundary into the uniformly distributed air ions. A 
step function compression in the magnetic field, narrower than the debris 
stream, is Initialized at the leading edge of the debris. The magnetic 
field is initially uniform on both sides of the compression. During the 
simulation runs, debris ions are continuously injected from the left 
boundary. This is accomplished by re-initializing the debris ion 
distribution in the left ghost cell each time step. 
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SAMSO/SZ 

POST OFFICE BOX 92960 
WORLDWAY POSTAL CENTER 
LOS ANGELES, CA 90009 
(SPACE DEFENSE SYSTEMS) 

OICY ATTN SZJ 

STRATEGIC AIR COMMAND/XPFS 
OFFUTT AFB, NB 68113 

OICY ATTN ADWATE MAJ BRUCE BAUER 
OICY ATTN HR T 

OICY ATTN DOK CHIEF SCIENTIST 

SAMSO/SK 
P.O. BOX 92960 
WORLDWAY POSTAL CENTER 
LOS ANGELES, CA 90009 

OICY ATTN SKA (SPACE COMM SYSTEMS) 
M. CLAVIN 

SAMSO/MN 

NORTON AFB, CA 92409 
(MINUTEMAN) 

OICY ATTN MNNL 

COMMANDER 

ROME AIR DEVELOPMENT CENTER, AFSC 
HANSCOM AFB, MA 01731 

OICY ATTN EEP A. LORENTZEN 

DEPARTMENT OF ENERGY 
LIBRARY ROOM G-042 
WASHINGTON, D.C. 20545 

OICY ATTN DOC CON FOR A. LABOWITZ 

DEPARTMENT OF ENERGY 
ALBUQUERQUE OPERATIONS OFFICE 
P.O. BOX 5400 
AL8UQUERQUE, NM 37115 

OICY ATTN DOC CON FOR D. SHERWOOD 


EG 8G, INC. 

LOS ALAMOS DIVISION 

P.O. BOX 809 

LOS ALAMOS, NM 85544 

OICY ATTN DOC CON FOR J. BREEDLOVE 

UNIVERSITY OF CALIFORNIA 
LAWRENCE LIVERMORE LABORATORY 
P.O. BOX 808 
LIVERMORE, CA 94550 

OICY ATTN DOC CON FOR TECH INFO 
DEPT 

OICY ATTN DOC CON FOR L-389 R. OTT 
OICY ATTN DOC CON FOR L-31 R. HAGER 

LOS ALAMOS NATIONAL LABORATORY 

P.O. BOX 1663 

LOS ALAMOS, NM 87545 


OICY 

ATTN 

DOC 

CON 

FOR 

J . 

WOLCOTT 

OICY 

ATTN 

DOC 

CON 

FOR 

R. 

F. TASCHEK 

OICY 

ATTN 

DOC 

CON 

FOR 

E. 

JONES 

OICY 

ATTN 

DOC 

CON 

FOR 

J. 

MALIK 

OICY 

ATTN 

DOC 

CON 

FOR 

R. 

JEFFRIES 

OICY 

ATTN 

DOC 

CON 

FOR 

J . 

ZINN 

OICY 

ATTN 

DOC 

CON 

FOR 

P. 

KEATON 

OICY 

ATTN 

DOC 

CON 

FOR 

D. 

WESTERVELT 

OICY 

ATTN 

D . 

SAPPENFIELD 



SANDIA LABORATORIES 
P.O. BOX 5800 
ALBUQUERQUE, NM 87115 

OICY ATTN DOC CON FOR W. BROWN 

OICY ATTN DOC CON FOR A. 

THORNBROUGH 

OICY ATTN DOC CON FOR T. WRIGHT 

OICY ATTN DOC CON FOR D. DAHLGREN 

OICY ATTN DOC COM FOR 3141 

OICY ATTN DOC CON FOR SPACE PROJECT 
DIV 

SANDIA LABORATORIES 
LIVERMORE LABORATORY 
P.O. BOX 969 
LIVERMORE, CA 94550 

OICY ATTN OOC CON FOR B. MURPHEY 

OICY ATTN DOC CON FOR T. COOK 

OFFICE OF MILITARY APPLICATION 
DEPARTMENT OF ENERGY 
WASHINGTON, D.C. 20545 

OICY ATTN DOC CON DR. YO SONG 
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OTHER GOVERNMENT 

INSTITUTE FOR TELECOM SCIENCES 
NATIONAL TELECOMMNICATIONS & INFO 
ADMIN 

BOULDER, CO 80303 


01 C Y 

ATTN 

A. 

JEAN (UNCLASS ONLY) 

01 CY 

ATTN 

W. 

UTLAUT 

01 CY 

ATTN 

D. 

CROMBIE 

01 CY 

ATTN 

L. 

BERRY 


NATIONAL OCEANIC & ATMOSPHERIC ADMIN 
ENVIRONMENTAL RESEARCH LABORATORIES 
DEPARTMENT OF COMMERCE 
BOULDER, CO 80302 
01CY ATTN R. GRUBB 
01CY ATTN AERONOMT LAB G. REID 


DEPARTMENT OF DEFENSE CONTRACTORS 
AEROSPACE CORPORATION 


P.O. BOX 92957 


LOS ANGELES, 

CA 

90009 

01 CY 

ATTN 

I. 

GARFUNKEL 

01CY 

ATTN 

T. 

SALMI 

01 CY 

ATTN 

V. 

JOSEPHSON 

01 CY 

ATTN 

S. 

BOWER 

01 CY 

ATTN 

D. 

OLSEN 


ANALYTICAL SYSTEMS ENGINEERING CORP 
5 OLD CONCORD ROAD 
BURLINGTON, MA 01803 

01CY ATTN RADIO SCIENCES 

AUSTIN RESEARCH ASSOC., INC. 

1901 RUTLAND DRIVE 
AUSTIN, TX 78758 

01CY ATTN L. SLOAN 
01CY ATTN R. THOMPSON 

BERKELEY RESEARCH ASSOCIATES, INC. 
P.O. BOX 983 
BERKELEY, CA 94701 

01CY ATTN J. WORKMAN 
01C Y ATTN C. PRETTIE 
01CY ATTN S. BRECHT 

BOEING COMPANY, THE 
P.O. BOX 3707 
SEATTLE, WA 98124 


01 CY 

ATTN 

G. 

KEISTER 

01 CY 

ATTN 

D. 

MURRAY 

01 CY 

ATTN 

G. 

HALL 

01 CY 

ATTN 

J. 

KENNEY 


CHARLES STARK DRAPER LABORATORY, INC. 
555 TECHNOLOGY SQUARE 
CAMBRIDGE, MA 02139 
01CY ATTN D.B. COX 
01CY ATTN J.P. GILMORE 

COMSAT LABORATORIES 
LINTHICUM ROAD 
CLAR<S8URG, MD 20734 
01C Y ATTN G. HYDE 

CORNELL UNIVERSITY 

DEPARTMENT OF ELECTRICAL ENGINEERING 
ITHACA, NY 14850 

01C Y ATTN D.T. FARLEY, JR. 

ELECTROSPACE SYSTEMS, INC. 

BOX 1359 

RICHARDSON, TX 75080 
01CY ATTN H. LOGSTON 

01C Y ATTN SECURITY (PAUL PHILLIPS) 

EOS TECHNOLOGIES, INC. 

606 Wilshire Blvd. 

Santa Monica, Calif 90401 
01C Y ATTN C.B. GAB8ARD 
01CY ATTN R. LELEVIER 

ESL, INC. 

495 JAVA DRIVE 
SUNNYVALE, CA 94086 
01CY ATTN J. ROBERTS 
01CY ATTN JAMES MARSHALL 

GENERAL ELECTRIC COMPANY 
SPACE DIVISION 
VALLEY FORGE SPACE CENTER 
GODDARD BLVD KING OF PRUSSIA 
P.O. BOX 8555 
PHILADELPHIA, PA 19101 
01 C Y ATTN M.H. BORTNER 
SPACE SCI LAB 

GENERAL ELECTRIC COMPANY 
P.O. BOX 1122 
SYRACUSE, NY 13201 

01 CY ATTN F. RE I BE R T 

GENERAL ELECTRIC TECH SERVICES 
CO., INC. 

HMES 

COURT STREET 
SYRACUSE, NY 13201 

01C Y ATTN G. MILLMAN 







GEOPHYSICAL INSTITUTE 
UNIVERSITY OF ALASKA 
FAIR8ANKS, AK 99701 

(ALL CLASS ATTN: SECURITY OFFICER! 

01CY ATTN T.N. DAVIS (UNCLASS ONLY! 
01CY ATTN TECHNICAL LIBRARY 
01C Y ATTN NEAL BROWN (UNCLASS ONLY) 

GTE SYLVANI.4, INC. 

ELECTRONICS SYSTEMS GRP-EASTERN DIV 
77 A STREET 
NEEDHAM, MA 02194 

01CY ATTN DICK STEINHOF 

HSS, INC. 

2 ALFRED CIRCLE 
3EDF0RD, MA 01730 

01CY ATTN DONALD HANSEN 

ILLINOIS, UNIVERSITY OF 
107 COBLE HALL 
150 DAVENPORT HOUSE 
CHAMPAIGN, IL 61820 

(ALL CORRES ATTN DAN MCCLELLAND) 

01CY ATTN K. YEH 

INSTITUTE FOR DEFENSE ANALYSES 
1801 NO. BEAUREGARD STREET 
ALEXANDRIA, VA 22311 
01C Y ATTN J.M. AEIN 
01C Y ATTN ERNEST BAUER 
01C Y ATTN HANS WOLFARD 
01C Y ATTN JOEL BENGSTON 

INTL TEL & TELEGRAPH CORPORATION 
500 WASHINGTON AVENUE 
NUTLEY, NJ 07110 

01CY ATTN TECHNICAL LIBRARY 

J AYCOR 

11011 TORREYANA ROAD 

P.O. BOX 85154 

SAN DIEGO, CA 92138 

01 CY ATTN J.L. SPERLING 

JOHNS HOPKINS UNIVERSITY 
APPLIED PHYSICS LABORATORY 
JOHNS HOPKINS ROAD 
LAUREL, MD 20810 

01 CY ATTN DOCUMENT LIBRARIAN 
01CY ATTN THOMAS POTEMRA 
01 C Y ATTN JOHN DASSOULAS 


KAMAN SCIENCES CORP 
P.O. ROX 7463 

COLORADO SPRINGS, CO 80933 
01CY ATTN T. MFAGHER 

KAMAN TEMPO-CEN T ER FOR ADVANCED 
STUDIES 

816 S Y A T E STREET (P.O DRAWER RG) 

SANTA BARBARA, CA 93107 
01CY ATTN DAS I AC 
01CY ATTN WARREN S. KNAPP 
01CY ATTN WILLIAM MCNAMARA 
01CY ATTN B. GAMBILL 

LINKABIT CORP 
10453 ROSELLE 
SAN DIEGO, CA 92121 

01CY ATTN IRWIN JACOBS 

LOCKHEED MISSILES 8 SPACE CO., INC 
P.O. BOX 504 
SUNNYVALE, CA 94088 
01CY ATTN DEPT 60-12 
01C Y ATTN D.R. CHURCHILL 

LOCKHEED MISSILES 8 SPACE CO., INC. 
3251 HANOVER STREET 
PALO ALTO, CA 94304 

01CY ATTN MARTIN WALT DEPT 52-12 
01CY ATTN W.L. IMHOF DEPT 52-12 
01CY ATTN RICHARD G. JOHNSON 
DEPT 52-12 

01CY ATTN J.B. CLADIS DEPT 52-12 

MARTIN MARIETTA CORP 
ORLANDO DIVISION 
'.U. BOX 5837 
ORLANDO, FL 32805 

01CY ATTN R. HEFFNEP 

M.I.T. LINCOLN LABORATORY 
P.O. BOX 73 
LEXINGTON, MA 02173 

01C Y ATTN DAVID M. TOWLE 
01C Y ATTN L. LOUGHLIN 
01CY ATTN D. CLARK 
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MCDONNEL DOUGLAS CORPORATION 
5301 80LSA AVENUE 
HUNTINGTON BEACH, CA 92647 
01CY ATTN N. HARRIS 
01CY ATTN J. MOULE 
01CY ATTN GEORGE MROZ 
01C Y ATTN W. OLSON 
01CY ATTN R.W. HALPRIN 
01C Y ATTN TECHNICAL 

LIBRARY SERVICES 

MISSION RESEARCH CORPORATION 
735 STATE STREET 
SANTA BARBARA, CA 93101 
01CY ATTN P. FISCHER 
01CY ATTN U.F. CREVIER 
01C Y ATTN STEVEN L. GUTSCHE 
01C Y ATTN R. BOGUSCH 
01CY ATTN R. HENDRICK 
01CY ATTN RALPH KILB 
01CY ATTN DAVE SOWLE 
01C Y ATTN F. FAJEN 
01 CY ATTN M. SCHEI8E 
01 CY ATTN CONRAD L. LONGMIRE 
01CY ATTN B. WHITE 
01CY ATTN R. STAGAT 

MISSION RESEARCH CORP. 

1720 RANDOLPH ROAD, S.E. 
ALBUQUERQUE, NEW MEXICO 87106 
01C Y R. STELLINGWERF 
01C Y M. ALME 
01C Y L. WRIGHT 


PENNSYLVANIA STATE UNIVERSITY 
IONOSPHERE RESEARCH LAS 
318 ELECTRICAL ENGINEERING EAST 
UNIVERSITY PARK, PA 16802 
(NO CLASS TO THIS ADDRESS) 

01C Y ATTN IONOSPHERIC RESEARCH 

PHOTOMETRICS, INC. 

4 ARROW DRIVE 
WOBURN, MA 01801 

01C Y ATTN IRVING L. KOFSKY 

PHYSICAL DYNAMICS, INC. 

P.O. BOX 3027 
BELLEVUE, WA 98009 

01C Y ATTN E.J. FREMOUW 

PHYSICAL DYNAMICS, INC. 

P.O. 3 0 X 10367 
OAKLAND, CA 94610 
ATTN A. THOMSON 


R R D ASSOCIATES 
P.O. BOX 9695 


MARINA 

DEL REY, 

CA 90291 

01 CY 

ATTN 

FORREST GILMORE 

01 CY 

ATTN 

WILLIAM B. WRIGHT 

01 CY 

ATTN 

WILLIAM J. KARZAS 

01 CY 

ATTN 

H. 

ORY 

01 CY 

ATTN 

C. 

MACDONALD 

01 CY 

ATTN 

R. 

TURCO 

0 1 C Y 

ATTN 

L. 

D e R A N D 

01 CY 

ATTN 

W. 

TSAI 


MITRE CORPORATION, THE 
P.O. BOX 208 
BEDFORD, MA 01730 

01C Y ATTN JOHN MORGANSTERN 
01CY ATTN G. HARDING 
01C Y ATTN C.E. CALLAHAN 

MITRE CORP 

WESTGATE RESEARCH PARK 
1820 DOLLY MADISON BLVD 
MCLEAN, VA 22101 
01CY ATTN W. HALL 
01CY ATTN W. FOSTER 

PACIFIC-SIERRA RESEARCH CORP 
12340 SANTA MONICA BLVD. 

LOS ANGELES, CA 90025 

01CY ATTN E.C. FIELD, JR. 


RAND CORPORATION, THE 
1700 MAIN STREET 
SANTA MONICA, CA <50406 
01CY ATTN CULLEN CRAIN 
01CY ATTN ED 3EDR0ZIAN 

RAYTHEON CO. 

528 BOSTON POST ROAD 
SUDBURY, MA 01776 

01C Y ATTN BARBARA ADAMS 

RIVERSIDE RESEARCH INSTITUTE 
330 WEST 42nd STREET 
NEW YORK, NY 10036 

01CY ATTN VINCE TRAPANI 








SCIENCE APPLICATIONS, INC. 

1150 PROSPECT PLAZA 
LA JOLLA, CA 92037 

01CY ATTN LEWIS M. LINSON 
01CY ATTN DANIEL A. HAMLIN 
01CY ATTN E. FRIEMAN 
C1CY ATTN E.A. STRAKER 
01CY ATTN CURTIS A. SMITH 

SCIENCE APPLICATIONS, INC 
1710 GOOD RIDGE DR. 

MCLEAN, VA 22102 
01CY J. COCKAYNE 
01C Y E. HYMAN 

SRI INTERNATIONAL 
333 RAVENSWOOD AVENUE 
MENLO PARK, CA 94025 
C1CY ATTN J. CASPER 
01C Y ATTN DONALD NEILSON 
01C Y ATTN ALAN BURNS 
01CY ATTN G. SMITH 
01C Y ATTN R. TSUNODA 
01CY ATTN DAVID A. JOHNSON 
01C Y ATTN WALTER G. CHESNUT 
01 CY ATTN CHARLES L. RINO 
01CY ATTN WALTER JAYE 
01CY ATTN J. VICKREY 
01 C Y ATTN RAY L. LEADABRAND 
01CY ATTN G. CARPENTER 
01 CY ATTN G. PRICE 
01CY ATTN R. LIVINGSTON 
01 CY ATTN V. GONZALES 
01CY ATTN D. MCDANIEL 

TECHNOLOGY INTERNATI ONAL CORP 
75 UIGGINS AVENUE 
BEDFORD, MA 01730 

01 CY ATTN W.P. BOQUIST 

TOYON RESEARCH CO. 

P.O. Box 6890 
SANTA BARBARA, CA 93111 
01 CY ATTN JOHN ISE, JR. 

01C Y ATTN JOEL GARBARINO 

TRU OFFENSE 8 SPACE SYS GROUP 

ONE SPACE PARK 

REDONDO BEACH, CA 90273 


01 CY 

ATTN 

R. 

K. PLEBUCH 

01 CY 

ATTN 

S. 

ALTSCHULER 

01CY 

ATTN 

D. 

DEE 

01 CY 

ATTN 

0/ 

STOCKWELL 


SNTF/1575 


VISIDYNE 

SOUTH BEDFORD STREET 
BURLINGTON, MASS 01803 
01CY ATTN W. REIDY 
01CY ATTN J. CARPENTER 
01CY ATTN C. HUMPHREY 

UNIVERSITY OF PITTSBURGH 
PITTSBURGH, PA 15213 

01CY ATTN: N. ZABUSKY 

DIRECTOR OF RESEARCH 
U.S. NAVAL ACADEMY 
ANNAPOLIS, MD 21402 
02CY 



















FILMED 


DTIC 






